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Preliminary results of a calculation of disconnected nucleon electromagnetic factors factors on the lattice are 
presented. The implementation of the numerical subtraction scheme is outlined. A comparison of results for 
electric and magnetic disconnected form factors on two lattice sizes with those of the Kentucky group is presented. 
Unlike previous results, the results found in this calculation are consistent with zero in these sectors. 



1. Introduction 

Fig.l illustrates a quark "disconnected dia- 
gram" in lattice QCD. These are present in al- 
most all hadron structure calculations and repre- 
sent one of the greatest calculational challenges 
in the field of lattice QCD. Experimental results 
from Thomas Jefferson National Laboratory are 
now closing on on these small effects in the form 
of the strange quark contribution to the nucleon 
electromagnetic form factors. 

The only previous published lattice results 
on electric and magnetic disconnected form fac- 
tors was in a 1998 paper by Dong, Liu and 
Williams H. They found that after chiral extrap- 
olation the disconnected parts gave 

G S M (0) = -0.36 ±0.20, 

\j\<r 2 s > n E \ w 0.25 < — > 0.4 fm. 

The last number is from monopole fits to the 
electric disconnected data. (Measured total is 

V\ < r2 >eL= °- 35 ± 0.03 fm). 
HAPPEX§ gives 

(G" E + 0.392G^)/(G A .///i p ) = 0.091 ± 0.54 ± 0.039, 

at q 2 = 0.477 (Gev/c) 2 . Similarly there is a 
less accurate number from SAMPLE^ at q 2 = 
0.1(Gcv/c) 2 : 

G^f (.l(Gcv/c) 2 ) = 0.61 ± 0.17 ± 0.21. 

Thus, experimentally there are no solid, separate 
numbers yet for disconnected G s M (q 2 ) or G s E (q 2 ). 

Theoretically, results for these quantities from 
different models are in poor agreement ||]. Pre- 
liminary results for disconnected electric and 




Figure 1. A generic nucleon disconnected dia- 
gram. The time axis is horizontal and the shaded 
circles represent nucleon sources or sinks. 

magnetic form factors will be presented here for 
the case where the disconnected quark has the 
same mass as the valence quarks. In order to ex- 
tract the strange quark contribution a chiral ex- 
trapolation of the different quark mass case must 
still be done. 

2. Noise and Subtraction Methods 

The basic method used to extract the signal is 
Z2 noise with perturbative noise subtraction. 
Consider Q such that 

< Tr{QX} >= 0, (1) 

where X is the noise matrix. One then has of 
course for general Q 

<Tr{{Q-Q)X} >=<Tr{QX}> . (2) 
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However, the variance will not be the same, 



V[Tr{(Q - Q)X}} ± V[Tr{QX}} 



(3) 



The natural choice for Q, which must mim- 
ick the off-diagonal components of Q to reduce 
the variance, is the perturbative expansion of the 
quark matrix [^-Q . This is given by 



M- X = I 



kP 



n 2 P 2 



k 3 P 3 



(4) 



where for Wilson fermions (I, J — {x, a, a}) 
(1 -J^UUx - a ll )S x , y+a J. 



(5) 



For a given operator, O, the matrix OM p 1 en- 
countered in < -0£>V> >^ -Tr^OMp 1 ) is not 
traceless. In other words, one must re-add the 
perturbative part, subtracted earlier, to get the 
full, unbiased answer. 

Local operators require perturbative correc- 
tions starting at 4th order in k (odd powers need 
no correction) and point-split ones require cor- 
rections starting at 3rd order (even powers need 
no correction). This procedure has been imple- 
mented exactly to 4th order (a single exact sub- 
traction) for the point-split conserved vector cur- 
rent on the lattice for the lowest 4 nonzero lattice 
momentums. 

3. Analysis Method 

A new analysis method was used for the dis- 
connected data analysis. At nonzero momentum, 
the basic quantity to consider is a ratio of nucleon 
three point and two point funxtionsQ], 

R G ( 2 ) M ) G p)(f, g r () 

where G( 3e,rra ) is the appropriate electric or mag- 
netic three point function and G^ 2 \t,q) is the 
nucleon two point function extending from the 
time origin to t with momentum q. The three 
point function comes from the correlation of the 
two point function and the loop data. Calling 
this Fourier transformed loop data L^ 6 '™) (q, t) = 
e~ lq ' x J(x, t), where J(x,t) is the appropriate 



component of the lattice self-contracted current, 
then generically 

G {3e ' m) {q,t',t) =< G (2e < m) (^0)L (e < m) (£',£) > 

- < G (2e < m) (i,0) >< L^ m \t',q) >, (7) 

where the average is over configurations. The sig- 
nal is always in the imaginary part of the lattice 
electromagnetic current J(x,t), which simplifies 
the analysis and improves error bars|g§[. 
The actual measurement used here is, 



t+i 

E 

f=i 



(R^ m \t,t',q)-R^ m \t-l,t',q)). (8) 



One can show for the appropriate components of 
G^ m \q,t' \t) and L^' m \t' ,q) that M^ rn \t - 



q) ->■ G E (q-), 



E+m 



G M (q ) for t > 1 in the 



electric or magnetic cases, respectively. The jack- 
knife technique was used to define correlated error 
bars and the time range t = 10 — 12, where the 
nucleon propagator becomes single exponential, 
was fit for all q. The measurement in Eq.(||) is 
represented graphically in Fig. 6. 

The use of the same time-summed loop back- 
ground in both terms in (8) leads to signifi- 
cantly smaller error bars. Eq.(|) implements the 
SESAM "plateau" strategy [10| by requiring the 
quark loop background to lengthen in time along 
with the nucleon propagators. Linear t fits of the 
individual terms in (8) nicely avoid systematic 
errors coming from "off-shell" or contact terms 
(t' ~ or t' ~ t, constants for t 3> 1 fll]]) since 
the loop time edges do not change their relation 
to the nucleon propagator time edges. However, 
the statistical errors are substantial. In contrast, 
the two terms in (8) could in principle possess dif- 
ferent off-shell contributions because the t' sum in 
the second term includes additional correlations 
with loop data two time steps beyond the nucleon 
sink (see Fig. 6). It has been explicitly verified 
that no discernable numerical correlations exist 
in such data, and so can be safely included. 

4. Results and Conclusions 

Two lattices sizes were used for comparison 
purposes: 20 3 x 32 and 16 3 x 24 at = 6.0. Re- 
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Figure 2. Disconnected nucleon electric form fac- 
tor at k — 0.152 on the 20 3 x 32 lattice as a 
function of q 2 in (GeV/c) 2 . 



Figure 4. Disconnected nucleon electric form fac- 
tor at ft = 0.152 on the 16 3 x 24 lattice as a 
function of q 2 in (GeV/c) 2 . Squares represent 
previous results from Ref. 0], circles are present 
results. Data here and in Fig. 5 separated in q 2 
for better visibility. 
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Figure 3. Same as Fig. 2 but for the disconnected 
nucleon magnetic form factor. 



Figure 5. Same as Fig. 4 but for the disconnected 
nucleon magnetic form factor. 
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Figure 6. Illustrating the two terms in the dis- 
connected amplitude measurement, Eq.(^). The 
time-summed quark loop data, symbolized by the 
underneath circles, is correlated with the appro- 
priate nucleon amplitude and multiplied by a ra- 
tio of two point functions. 



suits at k — 0.152 on both lattices will be exam- 
ined here (both use 100 configurations). 30 noises 
on the 20 3 x 32 lattice were used and 60 noises on 
the 16 3 x 24 lattices. As has been shown previ- 
ously || at fourth order subtraction for the elec- 
tromagnetic current, one gains a factor of about 
10 in the number of Z2 noises at essentially zero 
overhead. Thus these results are equivalent to 
about 300 unsubtracted noises on the large lat- 
tices and about 600 unsubtracted noises on the 
small lattices. 

Results are presented at k = 0.152 for the 
20 3 x 32 lattices in Figs. 2 and 3 above. At this k 
value there is no evidence for a signal in either the 
electric or magnetic sectors. In order to compare 
more directly with the results of Ref . jij , done on 
100 16 3 x 24 lattices with 300 unsubtracted noises, 
the calculation was repeated on identically sized 
lattices. Again, there is no evidence of a signal 
at nonzero momentum in the electric or magnetic 
sectors in Figs. 4 and 5. This finding is in dis- 
agreement with the Kentucky results ||||. The 



most serious disagreement is the electric sector, 
as can be seen in Fig. 4. 

The main message of these calculations to this 
point is that these correlations are extremely 
small QCD effects, both in the electric and mag- 
netic sectors. It will be crucial to find better 
methods which will extract the signals in a more 
efficient way. Isolating these signals and making 
a definitive prediction will be one of the great- 
est challenges facing hadronic lattice QCD in the 
years to come. 
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